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Abstract. We analyze the effects of the temperature on a bosonic Josephson junction 
realized with ultracold and dilute atoms in a double-well potential. Starting from 
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boson-boson interaction is suitably chosen thermal effects can increase the coherence 
visibility and reduce the on-site number fluctuation. 
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1. Introduction 

One-dimensional double-well potentials [1] confining ultracold and dilute bosons are 
the ideal arena to study the atomic counterpart of the Josephson effect [2]. Bosonic 
Josephson junctions (BJJs) have been widely explored at zero temperature. The lowest 
energetic state of a BJJ is reasonably described by that of the two-site Bose-Hubbard 
(BH) Hamiltonian [3]. This Hamiltonian sustains different ground-states depending on 
the coupling strength between the atoms. Then, varying the interatomic interaction 
makes possible to engineer macroscopic coherent states [U El El |8], macroscopic 
Schrodinger-cat states [HI HD1 EEH EE21 H31 E] and states characterized by high degrees 
of quantum correlations [HI [15]. The issue of which among the aforementioned states 
emerges in a BJJ is addressed by diagonalizing the underlying Hamiltonian and studying 
the (coherence) visibility of the interference fringes in the momentum distribution of 
the bosonic cloud, the quantum Fisher information, and the entropy entanglement as 
functions of the boson-boson interaction [T3j. It is worth to observe that the coherence 
visibility gives a measure of the importance of the single particle tunneling. On the 
other hand, the Fisher information is related in a very simple way to the on-site number 
fluctuation [16] which is the natural quantity to analyze the macroscopic quantum 
tunneling (MQT) of bosons across the barrier. Incidentally, notice that the single 
particle tunneling occurs on relatively short time scales, while the MQT takes place on 
times growing exponentially with the number of bosons. The situation studied within 
the two-site BH Hamiltonian framework is closely connected to that of MQT between 
the two classical Neel configurations of a two-site spin S easy axis antiferromagnet 
Hamiltonian [IT] which reasonably describes antiferromagnetic molecular rings, as Fe6 
and Cr 8 [181 CHj- In realistic situations, however, quantum systems work at finite 
temperature. For low dimensional geometries, effects due to thermal fluctuations are 
important [20]. Studying these effects on atoms in a double- well is very important as 
discussed experimentally in [21] and theoretically in [22] for spin-polarized fermions and 
in [23] for repulsive bosons. 

In this work we extend the studies presented in [13] and in [23]. In particular, we 
investigate the coherence and the on-site number fluctuation at finite temperature both 
in the attractive and in the repulsive regime by using as theoretical tool the two-site 
BH Hamiltonian. First we introduce the two-site BH model that we use to describe the 
Bose-Einstein condensate in a double-well. Then we diagonalize the BH Hamiltonian 
and study the coherence visibility and the on-site number fluctuation as functions of 
the interatomic coupling both at zero and at finite temperature. We point out that in 
the attractive regime the coherence visibility of the system exhibits an enhancement by 
increasing the temperature. We explain this effect by analyzing the coherence visibility 
of the thermally populated excited states. By using the Hellmann-Feynman theorem 
[23] and the time independent perturbation theory in the deep attractive regime we 
obtain an analytical formula for the coherence visibility pertaining to the two lowest 
BH Hamiltonian eigenstates. This formula is used to describe the afore mentioned 
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thermal enhancement. We perform a similar study for the on-site number fluctuation 
which exhibits a thermal softening in the attractive regime. Finally, we investigate the 
size effects on these quantities by varying the number of bosons in the system. 

2. The model Hamiltonian 

We consider N identical interacting bosons of mass m confined by a trapping potential 
Vtrapi?)- This potential can be realized by the superposition of an isotropic harmonic 
confinement in the transverse radial plane and a double- well potential (DWP) Vdw( x ) 
in the axial direction x. Then, V tr ap(j) is given by 

o 

7T? It) 

V trap (r) = V DW {x) + -£L + , (1) 

where u± is the trapping frequency in the radial plane. We assume that the transverse 
energy fkv± is much larger than the characteristic trapping energy along the x axis due 
to the potential Vdw{x). Then, the dynamics of the system is quasi one-dimensional 
(ID). 

In the following the system will be analyzed at finite temperature, T > 0. We 
assume that both k R T (k R is the Boltzmann's constant) and the characteristic boson- 
boson interaction energy are not much larger than the gap between the two states of the 
lowest doublet of the double- well linear problem [22], while they are much smaller than 
the gap between the first and the second doublet. In this case, only the two states in 
the lowest doublet will have non negligible occupancy. If the two wells are symmetric, 
the system will be described by the effective two-site Bose-Hubbard (BH) Hamiltonian 

® 

H = -J(a ] L a R + a R a L ) + ^{n L (h L - 1) + n R (h R - 1)) . (2) 

Here a k , a k (k = L, R, where L stays for left and R for right) are bosonic operators 
satisfying the usual commutation rules; h k = a k a k is the number of particles in the kth 
well; U is the boson-boson interaction amplitude, and J is the tunneling matrix element 
between the two wells. The total number operator N = Ul + n R commutes with the 
Hamiltonian (j2J). 

The spectrum of the Hamiltonian (J2j) is determined by solving the eigenvalues 
problem 

H\E i ) = E j \E i ) (3) 

for a fixed number N of bosons. Since the Hamiltonian H preserves the total number of 
particles, it can be represented by a (N + 1) x (N + 1) matrix in the Fock basis \i,N — i) 
(i = 0, N). In this ket, the left (right) index denotes the number of bosons in the left 
(right) well. For each eigenvalue Ej (j = 0, 1, ...,N) the associated eigenstate \Ef) will 
be of the form 

N 

m^^^N-i). (4) 

i=0 
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We assume that the coefficients c~p are real [14] . 

The macroscopic parameters of the Hamiltonian ([2]) are explicitly related to the 
atom-atom coupling constant g = 4:Trh 2 a s /m (with a s the s-wave scattering length), to 
the atomic mass m and the frequency u± of the harmonic trap (see, for example, |25j). 
The on-site interaction amplitude U is positive (negative) if a s is positive (negative), so 
that it may be changed at will by Feshbach resonance. Note that when a s < 0, to avoid 

the collapse, the system has to be prepared in such a way that the \U\ < — /t(w x wi)^ 2 , 
{25], with u x the trapping frequency of the single well. The hopping amplitude J is 
equal to (ei — eo)/2, where eo and t\ are the ground-state and the first excited state 
energies of a single boson in the double- well potential |22j . 

Note that at fixed number N of bosons, the ground-state of the system depends 
on the parameter £ = U/J [Sj. This parameter shall be used also throughout the 
present work to discuss the role of the boson-boson interaction in determining the 
system properties. In terms of ( the previous condition about the collapse thus reads 

|C| ^ j-j-jh(u x u)±) l l 2 . Then, fixed N and the radial frequency, the collapse can be 
avoided by suitably adjusting the parameters of the DWP. 



3. Analysis 

When the temperature is finite, the system is in a statistical mixture of states. We shall 
work in the canonical ensemble. In the basis of the eigenstates of the Hamiltonian (j2J) 
the mixed state density matrix p is 

N -I3E- 

p = Y. e —\ E ^\- ( 5 ) 

3=0 

Here Z is the partition function given by Z = J2f =0 e- pE *, where (3 = l/(k B T) with 
ks the constant of Boltzmann and T the absolute temperature. The thermal average of 
an operator A is then: 

(A) = J2—( E M\E ] ) , (6) 

3=0 

which, in the zero temperature limit, gets back the expectation value of A in the ground- 
state. 



3.1. Coherence visibility 

Recently, we have analyzed the coherence between the two wells of Vdw{ x ) a t zero 
temperature by calculating the coherence visibility [14J. At finite temperatures it is 
given by 

2 

a T = -v- T (a[a R ) , (7) 
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where the average (...) is the thermal average defined in Since we are assuming that 
the coefficients are real (see Sec. II), (o^or) = (o) R a L ) . The operator a} L a R (o^Sl) 
destroys a single boson in the right (left) well and creates it in the left (right) well. Then 
ay, see Eq. ([7]), characterizes the single particle tunneling through the barrier. The 
subsequent analysis will be then useful to understand how the temperature affects the 
tunneling of single bosons. 




Figure 1. (Color online). Coherence visibility a>T vs scaled interaction strength 
C = U/J at k B T = (solid line), k B T = 1.6J (dot-dashed line), and k B T = 2J 
(dashed line). Left panel: attractive bosons, £ < 0. Right panel (the horizontal axis is 
in logarithmic scale): repulsive bosons, ( > 0. Number of bosons N = 30. 

By using the expansion (jlj) in Eq. (j^J), ax reads 

9 N pE N 

«t = ^V^ * °& y/V + W-i) ■ (8) 

j=0 i=0 

Note that for i = N the squared root factor is zero. Therefore the summation in the 
right-hand side of Eq. (JSJ) is well defined even if c^ +1 is not defined in Eq. (j4j). 

In Fig. [1] we report ar as a function of the parameter £ = U/J both at zero 
and at finite temperature for N = 30 bosons. From the left panel, relative to the case 
of attractive bosons, we notice that, contrary to expectations, the coherence visibility 
at finite temperature is smaller than that at zero temperature only if the interaction 
strength is below a certain value. For interaction strengths larger than such a value the 
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coherence visibility exhibits an enhancement over its zero temperature value. We have 
numerically verified that this behaviour occurs when ksT > 0.5J. As we shall show, 
this corresponds to values of ksT of the same order of the gap between the two lowest 
eigenstates of the Hamiltonian (j2J). 

On the contrary, when the interactions are repulsive (right panel), we find that the 
coherence visibility at finite temperature (dot-dashed and dashed lines) is smaller or 
equal to its zero temperature counterpart (solid line) for any value of the interaction 
strength, in agreement with the experimental findings of Gati and co-workers |21j . 



0.02 
0.015 

0.01 
0.005 


-0.005 



1 


1 1 1 1 1 1 


1 1 1 1 1 1 1 1 




■ ■ ■ C= 0.1 


. <» ' ^ 




£=-0.2 


s' 

/' 






/ — 

/' 






/ 






/ 










/' 












/ 








y 




/ 






/ 


— ' 




/ 

/ 






/ 




1 


/ 

/ 




I 






1 


1,1,1, 


1,1,1,1,- 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 

k B T/J 

Figure 2. (Color online). Relative variation Sax of the coherence visibility vs scaled 
temperature ksT/J, calculated for scaled interaction strengths £ = —0.1 (dot-dashed 
line) and for £ = —0.2 (dashed line). Number of bosons N = 30. 



Due the unexpected behaviour of cit, we focus on the attractive regime and we 
investigate the relative change of the coherence visibility versus temperature for values 
of £ for which the coherence thermal enhancement is expected to occur. In particular, 
we study the relative variation 

CtT — «o 



5ai 



(9) 



with a the coherence visibility at T = and report the results in Fig. [2j When 
C = —0.1 (dot-dashed line), 5ax is negative up to a given temperature, after that it is 
positive. For ( = —0.2 (dashed line), the relative change of the coherence visibility is 
always positive when the thermal energy is sufficiently high. The positive values of 5oit 
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for ( = —0.1 (dot-dashed line of the Fig. [2]) are larger than those at £ = —0.2 (dashed 
line of Fig. |2J). 
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Figure 3. Relative variation Scut of the coherence visibility vs number of bosons N, 
calculated at k B T = 2 J for ( = -0.1. 

The thermal enhancement of the coherence visibility can be related to the fact that 
at T > bosons may populate excited states of the Hamiltonian (j2J) which have a larger 
coherence than the ground-state. To support this interpretation we use a simplified 
model in which we suppose that fc^T ~ (Ei — Eq) so that only the two lowest eigenstates 
of the Hamiltonian fl2]) are populated. At this stage, the thermal coherence visibility 
may be written as 



«i e 



Here (Xj 



l + e-i'l + e-T" < 10 > 
(Ej\a L api\Ej) is the visibility of the jth excited state of the two-site BH 
Hamiltonian. We, then, have to prove that the right-hand side of Eq. ( JTUl) is larger 
than a , i.e. 

a 1 — a > . (11) 

We evaluate ot\ — «o by jointly exploiting the Hellmann-Feynman (HF) theorem [21] 

and the time independent perturbation theory when U/J — > — oo which provides 

results in excellent agreement with those deriving from the exact diagonalization also 

for moderately small values of \(\ as found in [T3] for N = 30. As observed in [H], the 

HF theorem allows to write: 

ldE j 



ft; 



N dJ 



;i2) 
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so that: 

l,dE dE, 

ai - ao = iv ( 9J- 97 } - (13) 

In the limit U/J — > — oo, the hopping operator in the Hamiltonian ([2]) can be treated 
within the perturbation theory. We calculate E and Ex at the first non vanishing order, 
and get for ot\ — «o the following expression: 

4(iV+l) 

ai ~ ao = -N(N-1)(N-3)C ' ( 4) 

which is positive when N > 3 when the interaction is attractive. 

We can explicitly estimate the temperature at which thermal effects affect the 
coherence visibility. This will happen when k R T ~ (E\ — Eq). From perturbation 
theory one has: 

Equating this value to k R T, in the limit of N ^> 1 we obtain: 

J ^(UNl + UNk s T (16) 

At temperature T, the interaction energy —UN/2 ~ ksT, so that we get 

J — ksT . (17) 

It is fair to remark that, when present, the increase of the coherence visibility over 
its zero temperature value is quite small and it depends on the number of bosons N. 
This is shown in Fig. [31 where we have plotted 5 ax with olt evaluated at k R T = 2 J as a 
function of N for £ = —0.1. The vanishing of 5ax with ax when N increases is confirmed 
also for higher temperatures and larger attractive strengths. From Fig. |3l moreover, 
we can see that the maximum gain in visibility (~ 0.04) is reached when N = 34. Such 
a moderate iV corresponds to a rather significant shot noise and experimentally it is 
necessary to gain sufficient statistics to demonstrate the contrast enhancement. 

Summarizing, we have found that for attractive bosons when both the temperature 
of the system and the boson-boson interaction become sufficiently large, the coherence 
visibility exhibits a thermal enhancement which, however, vanishes for large N. Since 
the system under investigation is a closed system, its coherence properties are only 
determined by the interplay between the tunneling coupling energy (resulting from the 
overlap of the wave functions localized in the two wells) and the localization energy due 
to the the interaction between the particles. This has been shown in [5] for the case 
of repulsive interactions and can be seen also for attractive bosons by evaluating the 
expectation value E of the BH Hamiltonian (j2J) in the state [H] 

\QC) = \CS) L ®\CS) R , (18) 

where \CS)k (k = L,R) is the coherent state [27] describing the Bose-Einstein 
condensate in the kth well. Thus, the expectation value E (up to a constant) reads 

n 



N 2 U 



E = (QC\H\QC) = -N Jy/1 - z 2 cos </> + —— z 2 , (19) 
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where z = (Nl — Nr)/N is the fractional population imbalance between the two wells 
and <f) = Or — 9l is the phase difference between the two condensates. Under the 
hypothesis of sufficiently small population imbalance, the energy (1191) becomes 

E = -Ej cos <p + E c (An) 2 , (20) 

where An = (N L — N R )/2, Ej = NJ is the Josephson coupling energy, and E c = U 
is the localization energy due to the boson-boson interaction. When both temperature 
and boson-boson attraction are sufficiently large, Ej overcomes the localization energy 
E c by producing a visibility increasing with the temperature. 



3.2. On-site number fluctuation 

In [13] we have also analyzed the quantum Fisher information (QFI) at zero temperature. 
The QFI is used as a probe of the emergence of the NOON state - written as a symmetric 
combination of the Fock states \N, 0) and |0, N) - and is related to the on-site number 



fluctuation [T6] . 

In the following we shall investigate the thermal average of the squared of the on-site 
number fluctuation (Ahk)% (k = L,R as usual). According to Eq. (jHJ): 



N e~^ 



(An k ) 2 T = ^——(EjlnllEj) 



j=0 

N -f3E- 2 

E^^i^i^)) • ( 21 ) 

j=0 

The operator = a\hk counts the number of bosons in the kth well. Given a state of 
the bosonic junction, the squared variance (|2"T|) measures the deviation of the kth well 
population from its expected value, such a deviation being due to the particles tunneling 
between the two wells. Then, (An^jy can be used to characterize the collective transfer 
of bosons across the central barrier of Vdw(%)- in the following we investigate the 
thermal effects on this macroscopic quantum tunneling (MQT). 

To fix the ideas, let us focus on the number fluctuation in the left well, k = L. In 
terms of the coefficients of the expansion (j3J), (An L )x = /t is given by 

N P -PE 3 N 

j=0 i=0 

N _ BE . N „ 

(22) 

We have studied fr as a function of the scaled interaction parameter ( — Uj J by 
varying the temperature, see Fig. HI The solid line of this figure represents the on-site 
number fluctuation at T = 0, the dot-dashed and the dashed lines give the same quantity 
calculated, respectively, at k B T = 1.6J and k B T = 2J. From the plots of Fig. HJ we 
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Figure 4. (Color online). On-site number fluctuation fx for the left well population 
vs scaled interaction £ = U/J at fc^T = (solid line), fcgT = 1.6 J (dot-dashed line), 
and fcgT = 2 J (dashed line). Left panel: attractive bosons, £ < 0. Right panel (the 
horizontal axis is in logarithmic scale): repulsive bosons, £ > 0. Number of bosons 
N = 30. 



see that at finite temperature fx is enhanced over its zero temperature value for both 
attractive (for sufficiently weak boson-boson interactions) - left panel - and repulsive 
bosons, right panel. When the bosons are attractively interacting, it exists a value of 

say | £|, above which the intra- well number fluctuations at T > are smaller than 
those at T = 0. The strength of this interaction is very close to that for which the ff of 
the ground-state has a vanishing second derivative with respect to (. It is interesting to 
observe that when T = and the interaction strength approaches \(\, the system begins 
to lose its coherence and the junction evolves towards the self-trapping regime |14j . 

We now analyze the relative change of fx both in the attractive and in the repulsive 
regime. To this end, we study 5 fx defined by: 

(An L ) 2 T - (An L ) 2 

5fT ~ (A^l (23) 

as a function of the temperature in correspondence to different interaction strengths. 
We have carried out this analysis for N = 30 and N = 100, see Figs. [5] and El 

On the repulsive side (top panels of Figs. |5]and[6]), the relative change (123]) increases 
with the temperature. The influence of the interatomic interaction is to reduce such a 




Figure 5. (Color online). Relative variation 5fr of number fluctuation vs scaled 
temperature kgT/ J. Number of bosons N = 30. 



growth with the increasing of £. The middle and bottom panels of Figs. and M are 
obtained for attractive bosons. From the middle panels (see the dashed lines therein) 
we observe that when the boson-boson attraction becomes sufficiently strong, a change 
in the concavity of 5 fx takes place. By looking at the bottom panels we see that a 
further increasing of the interatomic attraction produces a maximum and a minimum 
(C ~ -0.087 for iV = 30 and ( ~ -0.023 for N = 100) in 5f T that disappear for 
sufficiently high values of £. At this point (£ ~ —0.12 for N = 30 and £ ~ —0.025 for 
= 100) Sfr is zero up to a given temperature, after that it decreases against ksT/J, 
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Figure 6. (Color online). Relative variation 5fr of number fluctuation vs scaled 
temperature fcsT/J. Number of bosons N = 100. 



as shown by the solid lines of bottom panels of Figs. and El From these solid lines, it 
can be observed that the greater is N the weaker is the attraction at which the MQT 
thermal softening sets in. 

To explain the thermal softening of fx we follow the same path of reasoning used for 
the coherence visibility in the previous subsection . We suppose that ksT ~ (Ei—E ) so 
that only the two lowest eigenstates of the Hamiltonian ((2]) are populated. The thermal 
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on-site number fluctuation fa may be written as 

fo fx e- 1 



h 



l + e- 1 1 + e" 1 

\2| 



(24) 



where fj = (Ej\(Arik) \Ef). We, then, have to find the conditions under which the 
right-hand side of Eq. 0240 is smaller than / , that is 



fi - fo < . (25) 

To calculate fi — fo we use, again, the time independent perturbation theory and the 
HF theorem to get [H 



fi 



and therefore: 



fi-f 



dU ' 2 x ~ 2 
d{E 1 - E ) 



dU 



(26) 



(27) 



In the deep attractive regime, we can treat the hopping operator within the perturbation 
theory and calculate Eq and E\ to the first non vanishing order to get f\ — f§. 

2(iV + l) 



A-/o = -( 



N+ 1 + 



(N- l)(iV-3)C 2 / ' 



(28) 



which is negative for N > 3. 
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Figure 7. (Color online). Relative variation 6 fx, of the number fluctuation, vs number 
of bosons TV, calculated at fcsT = 2 J. Dashed line: £ = —0.01. Solid line: £ = 0.1. 
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Figure 8. Relative variation 8 fx, of the number fluctuation, vs number of bosons TV, 
calculated at k B T = 2J for ( = -0.1. 

We have also analyzed the influence of the number of bosons both on the thermal 
enhancement and on the thermal softening of the on-site number fluctuation. Let us 
start with the case 5 fx > shown in Fig. It can be pointed out that when the 
bosons interact attractively (dashed line) the relative variation of fa increases with N 
by reaching a maximum, after that it decreases at least up to the number of bosons 
that we have considered. On the other hand, for repulsive bosons (dotted line), 8 fx 
decreases by increasing N; in correspondence to the same number of particles, the 
thermal enhancement of fx is more important for attractive bosons. The case 5 fa < 
is reported in Fig. [HJ It can be observed that 6fy increases by increasing the number 
of particles in the system. 

Let us, now, write down some conclusive remarks. The new predicted effects - 
coherence thermal enhancement and on-site fluctuation number thermal softening - 
would manifest at temperatures of about J/ks (see Fig. 2 and Fig. 6). In the 
experiments performed to detect the scaling behavior of the coherence factor in a bosonic 
Josephson junction, the temperature varies in the range: 50nK < T < 80nK, see, for 
example, [21]. It is fair to note that to observe the two mentioned phenomena in the 
above temperatures range, very large values of the hopping amplitude (corresponding 
to Vdw with very low central barriers) - 6.6kHz < J/h < 10.5kHz - would be required. 

Possible limitations to the feasibility of the scenario analyzed in the present work 
are related to the assumption of the thermal equilibrium that we have made, see Eq.([5]). 
Only a sufficiently weak coupling between the confined bosons and an external thermal 
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bath does not affect the physical properties of the system. For this reason, we expect that 
the coherence thermal enhancement and on-site fluctuation number thermal softening 
may be observed in mesoscopic systems only if they are very weakly coupled to the 
environment. When this condition is not met, the effect of the number fluctuation has 
to be considered, in particular for a moderately small number of bosons. This problem 
is currently under study. 

4. Conclusions 

We have considered ultracold and dilute bosons in a one-dimensional double-well poten- 
tial. By considering such a system at finite temperature, we have studied the thermal 
effects both on the single particle and on the macroscopic quantum tunneling. We have 
carried out this analysis by studying the coherence visibility and the on-site number 
fluctuation as functions of the temperature for different interaction strengths. We have 
pointed out that the thermal effects can increase the coherence visibility and reduce the 
on-site number fluctuation when the interatomic interaction is suitably tuned. We have 
explained the coherence thermal enhancement by analyzing the coherence visibility of 
the thermally populated excited states. By employing the Hellmann-Feynman theorem 
and the time independent perturbation theory, we have explicitly evaluated the coher- 
ence visibility pertaining to the two lowest eigenstates of the two-site Bose-Hubbard 
Hamiltonian. A similar approach was been followed to calculate the on-site number 
fluctuations of the lowest states of the two- mode Hamiltonian and justify the observed 
thermal softening of the intra-well fluctuation. We have investigated size effects on the 
two mentioned quantities by varying the number of bosons in the system. 
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